Hierarchical energy relaxation in mesoscopic tunnel junctions: 
Effect of nonequilibrium environment on low-temperature transport 
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We develop a theory of far from the equilibrium transport in arrays of tunnel junctions. We find 
that if the rate of the electron-electron interactions exceeds the rate of the electron-phonon energy 
exchange, the energy relaxation ensuring the charge transfer may occur sequentially. In particular, 
cotunneling transport in arrays of junctions is dominated by the relaxation via the intermediate 
bosonic environment, the electron-hole excitations, rather than by the electron-phonon mechanism. 
The current- voltage characteristics are highly sensitive to the spectrum of the environmental modes 
and to the applied bias, which sets the lower bound for the effective temperature. We demonstrate 
that the energy gap in the electron-hole spectrum which opens below some critical temperature T* 
due to long-range Coulomb interactions gives rise to the suppression of the tunneling current. 

PACS numbers: 72.10.-d, 73.23.-b, 73.63.-b, 74.50. +r 



Transport in mesoscopic tunnel junctions is ensured by 
the energy exchange between the tunneling charge car- 
riers and energy reservoirs: since the electronic energy 
levels at the banks of the mesoscopic junctions are, in 
general, different, the tunneling is impossible unless there 
is subsystem of excitations capable of accommodation of 
this energy difference @, §,§, § |, g, @j. Intense stud- 
ies of nano-structured and disordered systems including 
Joscphson junctions mesoscopic superconductors 
patterned superconducting films |l(J, highly disordered 
superconducting and semiconducting films jll], |l2|, [l3|, Q 
reveal a prime importance of the out- of- equilibrium prop- 
erties of an environment to which the tunneling charge 
carriers relax the energy. Notably, the relaxation pro- 
cesses can be mediated not only by phonons but by 
the energy exchange with the electromagnetic environ- 
ment |],Ti |, [I |, @, H 01 and with the electron-hole 
(c-h) pairs generated by the tunneling carriers ji~7|, [l8|| . 
The energy relaxation in mesoscopic tunnel junctions in 
the case where the energy exchange between the tun- 
neling carriers and the electromagnetic and/or electron- 
hole reservoir, l/r cnv _ c , is comparable to the rate of the 
energy loss to the phonon thermostat, l/T env -»bath) was 
analyzed in jy^. In this Letter we develop a general 
approach to the description of the strongly nonequilib- 
rium processes where l/r cnv _ c 3> l/V e nv->bath and show 
that the energy relaxation enabling the tunneling cur- 
rent occurs in two stages: (i) The energy relaxation from 
the tunneling charges to the intermediate bosonic modes 
(electromagnetic or electron-hole excitations) which we 
hereafter call the environment] and (ii) The energy trans- 
fer from the environment to the phonon thermostat, to 
which we will be further referring as to a bath. 

We demonstrate that the transport is controlled by the 
first stage and is thus critically sensitive to the spectrum 
of the environmental modes. At the same time, the pass- 



ing current drives the environment out of the equilibrium, 
and the environment spectrum and effective temperature 
may become bias-dependent themselves. We derive the 
coupled kinetics equations for the charge carriers and out- 
of-equilibrium bosonic environment and apply our tech- 
nique to tunneling transport in large arrays of normal 
and superconducting junctions. 

A single junction - First, we consider a tunnel junc- 
tion between two bulk metallic electrodes biased by the 
external voltage V, see Fig. Ek. A general formula for the 
tunneling current reads: 



I = e 



(r-r) 



(1) 



where T ( T ) is the tunneling rate from the left (right) 
to the right (left), and, for a single junction, 
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/W(l-/<? ) )P<(e-e'), 
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where f^ 1,2 ^ are the electronic distribution functions 
within the electrodes, P < (e) is the probability for the 
charge carrier to lose the energy E to the environment, 
and R T is the bare tunnel resistance. The backward 
scattering rate, V oc / ££ , / e (2) (l - f i e P)P < (e - e')- If 
an intermediate environment is absent and the relax- 
ation is provided by the phonon bath, then P < (e) = 5(e) 
and Eq. (||) reproduces the conventional Ohm law. The 
quasi-equilibrium situation where the distribution func- 
tions of the environmental modes N u are Bose distribu- 
tions parameterized by the equilibrium temperature was 
20 1 . In a general, far from the equilibrium 
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FIG. 1: (Color online.) (a) The effective circuit for the 
tunnel junction subject to bias V and with the environ- 
ment having the impedance Z. (b)-(c) Diagrammatic expan- 
sion of P < to the first and the second orders in p respec- 
tively. The solid lines represent propagation of electrons, the 
dashed lines denote the environment excitations. The ver- 
tex with the two electron lines and one dashed line carries 
GtP{oj)/uj factor, the two dashed- lines vertex corresponds to 
G T p(w)p(w') /(loco') . 



case, we find: 
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F(ui) = [N^ + (1 + N^e 



(3) 

(4) 
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The nonequilibrium distribution function N u is defined 
by the kinetic equation with the scattering integral 
describing the energy exchange between environmental 
modes and tunneling electrons. Terms proportional to 
the Nu and 1 + N u correspond to the absorbed and emit- 
ted environmental excitations respectively. The combina- 
tion B u = 1 + 2N UJ is the kernel of the time-independent 
contribution to J describing the elastic interaction of the 
tunneling electron with the environmental modes and 
having the structure of the Debye- Waller factor. In an 
equilibrium, N u reduces to the Bose-function and the 
functional P < recovers the result by Rcf. p(j|] . The 
spectral probability of the electron-(electromagnetic) en- 
vironment interaction is p(u>) = Ke[Z t (u>)]/ R Q , where 
Z t is the total circuit impedance and R Q is the quan- 
tum resistance ||. Proceeding analogously to Ref. [pTJ , 
one finds the spectral probability corresponding to the 
electron-environment interaction within the each elec- 
trode as pi(io) = 21m J Ui/(Diq 2 - i^) 2 , i = 1,2, and 
that for the interaction across the junction, /9i 2 (cj) = 
— 21m J Ui2/[(Diq 2 — iu;)(.D 2 (? 2 — ioS)], where -Di( 2 ) are 
diffusion coefficients within respective electrodes, and 
U\{2) are the dynamically screened Coulomb interactions 
within (across) the electrodes. The form of p(u>) de- 
pends on the structure of the environmental excitations 
spectrum and, thus, on the external bias. The latter is 
especially important in the array of highly transparent 
junctions where p{uj) is different for elastic and inelastic 
processes jl], p2|]. In particular, for the e-h environment 



with the constant U, one should cut off the integral at 
q = \jT c ff/D, when calculating p(0), where T c s is the 
(bias dependent) effective temperature of the environ- 
ment which we determine below. This allows us to for- 
mulate a recipe: if in an equilibrium p = p(lo, T) then in 
an out-of-equilibrium state p = p(u>,T e g). 

To close the set of formulas (0)-(||) one has to add 
the kinetic equations (KE) for the boson distribution 
functions N u . To derive these KE we use a semi- 
phenomenological kinetic approach of |23[ ] and express 
the current of Eq. (|l|) through the electronic distribution 
function as I = J vi[dfel /dt]. Here df^ /dt = I co i, 
where I C oi is the collision integral describing the evolu- 
tion of the electronic distribution function due to en- 
ergy and/or momentum transfer processes. Expanding 
further P < with respect to p we obtain, in the zero 
order in N u , the collision integral in a form ij^ — 

-JWMf^il - /i 2 2) ) - /i 2 2) (l - f^Md - e 2 > 2 de 2 , 
where W12 = l/ui^Rr is proportional to the bare prob- 
ability for an electron to be transmitted from one lead to 
the other. In the first order 

dfei^ f ( P \ 

— ^— = - diov u) v 2 de2 W12X 

dt J \ luv,,, / 



5(^12 - + !)(!- fg>) (1 - f£>)N u fg>]+ 



(2)^ 



P(2)l 



S(e 12 + w)\f^N u (l /4 2 )) - (1 - fg))(N u + l)/ £ (2) ] 



where ei 2 = e\ — e 2 and z/ w is the density of environmental 
states 0] . The structure of I^J is identical to that of the 
electron-phonon scattering integral in metals |]23f , where 
Nu would stand for the phonon distribution functions. 
The quantity p/u)u u is proportional to the probability of 
the electron-environment scattering. 

The collision integral dual to and describing the 
evolution of N u is derived analogously, and the resulting 
kinetic equation is: 



dNu 
dt 



VuiR-T 



[Nuil + nu,) - (l + JVuH 



(6) 



where A is the numerical factor of order of unity, n w is 
the electron-hole pairs distribution function. The scat- 
tering integral in Eq.(|^) is also identical by its struc- 
ture to the phonon-electron scattering integral in met- 
als [^3) . For the electron- hole environment (i = 1,2 
label the electrode in which the pair is located), one 
has riu = (1/ui) f ,fil (1 — this agrees with the 

results of Ref. fl25| where the nonequilibrium boson dis- 
tribution function is equivalent to our 1 + 2n u . If elec- 
trons and holes belong to different electrodes, n w = 
{2co)- 1 ff^af j (l - d* being the Pauli-matrix. 

From ra) one estimates the rate of the energy exchange 
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between the environment and the tunneling electrons 
as: l/r cnv _ c = p(uj)/u u R T . Now one has to compare 
1 /T cnv _ c with the rate of the interaction of the environ- 
ment modes with the (phonon) bath, 1 / T env ^ a th(u>) . For 
the electron-hole environment, l/T e nv-»bath(w) is deter- 
mined from Eq. (|^) to which the electron-phonon scatter- 
ing integral is added. If T onv ^bath ^ 7"cnv-c, the two- 
stage relaxation takes place and the characteristic en- 
ergy transfer from tunneling current is to ~ max{T c , V}, 
where T e is the electronic temperature in the leads. The 
electromagnetic environment mediates the two stage re- 
laxation in the case where Ohmic losses occur in a LC 
superconducting line and are small S]. To take a typ- 
ical example, in aluminum mesoscopic samples r en v-e — 
10 _8 sec and r env _,bath = 10~ 6 scc jm], so the conditions 
for the two-stage relaxation are realized. In this case the 
distribution functions, N u , of the environmental modes 
become nonequilibrium and are determined from the con- 
dition that the collision integral of the environmental 
modes with the e-h pairs accompanying the current flow 
becomes zero, then Eq.(|^) yields N u = n^} 2 \ If T <C V, 
then Nu can be approximated by the Bosc-function with 
some effective temperature T e g at u> < V = T e ff and 
N u = at uj > T e g (the emission of the excitations with 
the energy larger than V is forbidden), and 



T eS = lim luN^ 



0.5Fcoth(y/2T) 



(7) 



Thus the system with the environment well isolated from 
the bath cannot be cooled below T e g. 

Equations ([!])- @ give the full description of the ki- 
netics of the tunneling junction in a nonequilibrium en- 
vironment. To derive the I-V characteristics we find 
N u = n£ 12 ^ and plug it into Eqs.(Q)-(Q). Introducing 
the parameters g = p(0) and A, the characteristic fre- 
quency of the p(lu) decay [for the Ohmic model p0[ , 
p = g^ 1 /{l + (uj/K) 2 } and A/g is of the order of the 
charging energy of the tunnel junction], we find: 



Rr V 



(8) 



in the interval T < V < A, where T eS ~ V. Note 
that /(V) given by Eq.( p|) differs from the power law 
dependence obtained in |20|] for T c = T ff = 0. This 
shows that tuning the environment one can manipulate 
by the tunnel junction I(V) (the gating effect). At high 
voltages, V 3> A, one finds 



I{V) ~ (V - Aoo)/i? T , 



(9) 



A, = »J'(0) = 2/ X3 ^[l + Ad out) - M in) ] * 
A^ ) ln(A/min{T e ,T cnv }), where = A^N^ = 

N^} ~ A/g, since at V > A, A^ out) ~ A/w > A# n) . 
Arrays of tunnel junctions- Extending Eq. (|^) onto an 
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FIG. 2: (Color online.) (a) The single electron two-islands cir- 
cuit. (b)-(e) Diagrams describing the forward inelastic cotun- 
neling rate. The "up" arrows stand for the e-h pairs excited 
during the cotunneling and the "down" arrows correspond to 
the recombination of the e-h pairs. The vertices shown by 
boxes are proportional to the probability of an elemental e-h 
pair excitation, p(u>)/cj. 



array comprised of junctions one finds 
/ N 



Rk 



Att 2 R, 



S 2 I dede'h{e)[l-h{e')]P{e-e'), 

(10) 



where 



P(E) 



dtexp(iEt) 



n [ 



Ar(in) iwt 



j p{u) 
alo x 



+ (1 + Nl°f)e 



u,3 



j<N-l 



Here S 



-(JV-i) atN 



N N /{N- 1)!, and E c 



(11) 



2 /2C is 



the Coulomb charging energy of a single junction (C is 
a single junction capacitance) and for the Cooper pair 
transport e — > 2e. Eqs. (|l0|) , (pd|) were derived in a first 
order in tunneling Hamiltonian. Shown in Figj2| is a dia- 
grammatic representation of Eq.(]ll|) for N = 3. 

A generalization of the results obtained for a single 
junction including the structure of the collision integral 
and the concept of the effective temperature Eq. (^) , onto 
large arrays is straightforward. As long as tempera- 
tures are not extremely low p7[ , the charge transfer in 
large arrays is dominated by the inelastic cotunneling 
and the two-stage energy relaxation. The tunneling car- 
riers generate e-h pairs |l7], 18 serving as an environ- 
ment exchanging the energy with the tunneling current 
and then slowly losing it to the bath. It is instructive 
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to consider a two-dimensional array of superconducting 
tunnel junctions. On the distances L < A = y/ C/Cq, 
where Co is the capacitance of a single junction to the 
ground, the Coulomb interaction between charges is log- 
arithmic. If the size of an array does not exceed A, 
the e-h plasma comprising the environment experiences 
the charge Berezinskii-Kosterlitz-Thouless (BKT) tran- 
sition [§6[ |7) at T = T BKT ~ E c |27|. This implies 
that at T ~ E c the energy gap T* opens in the spec- 
trum of unbound electrons and holes and, as a result, 
p(uj) vanishes in the interval < uj < T*. One than 
sees from Eq. ( |lCi| ) that opening the gap suppresses both 
Cooper pairs- and normal quasiparticle currents in the 
superconducting tunneling array at T < E c . Analyzing 
contribution from higher orders into cotunncling process, 
one finds that the current suppression holds in all orders. 
This picture applies to the films close to superconductor- 
insulator transition (SIT) Q. Indeed, near the SIT the 
dielectric constant e of the film diverges 1 28 and on the 
distances L < ed, where d is the film thickness, the 2D 
e-h environment experiences the BKT transition. Thus 
opening the gap in the electron-hole spectrum due to 
long range Coulomb effects and the resulting suppression 
of the tunneling current offers a microscopic mechanism 
for the insulator-to-superinsulator transition [^9| [30| . 

Two notion of the two-stage relaxation is a key to 
resolving the controversy of the variable range hop- 
ping (VRH) conductivity in both doped semiconduc- 
tors |n| |2| and disordered superconducting films [ |L2| : 
the observed universal pre-exponential factor indicates 
that the energy relaxation is due to electron-electron (e- 
e) rather than the electron- phonon interactions. On the 
other hand, according to |33L Q e-e relaxation cannot en- 
sure a finite conductivity below the so called many-body 
localization temperature J^|. The sequential relaxation 
of hopping electrons via the e-h environment, which fur- 
ther transfers energy to the phonon bath implies that the 
prefactor in hopping conductivity is indeed proportional 
to e 2 /h. 

In conclusion, we have developed a quantitative de- 
scription of the highly noncquilibrium tunneling trans- 
port in arrays of tunnel junctions in the limit l/r cnv _ c 3> 
1 /^env^bath and demonstrated that the low-temperature 
relaxation ensuring the tunneling current occurs via an 
intermediate electromagnetic and/or e-h pairs environ- 
ment. We argued that the onset of the gap in the spec- 
trum of environmental excitations suppresses tunneling 
current. In particular, the gap due to Coulomb interac- 
tions in superconducting arrays can offer a microscopic 
mechanism for the insulator-superinsulator transition. 
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